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Abstract. The possibility is investigated that the Ha luminosity at the nebular epoch may be an additional 
indicator of ^''Ni mass in type II supernovae with plateau (SNe IIP), on the basis of available photometry and 
spectra. We first derive the ^''Ni mass from the My magnitude on the radioactive tail using a standard approach. 
A confirmation of the correlation between ^''Ni mass and plateau My magnitude found recently by Hamuy (2003) 
is evident. There is strong evidence of a correlation between steepness of the V light curve slope at the inflection 
time and the ^^Ni mass. If confirmed, this relation may provide distance and extinction independent estimates 
of the amount of ^®Ni in SNe IIP. We then apply upgraded radioactive models of Ha luminosity at the nebular 
epoch and claim that it may be a good indicator of ^''Ni, if mass, energy and mixing properties vary moderately 
(within factor ~ 1.4) among SNe IIP . This method of the '^'^Ni mass determination from Ha luminosities yields 
results which are consistent with the photometric mass of ''^Ni mass to within 20%. This result also implies that 
the parameters of SNe IIP events (mass, energy and mixing properties) are rather similar among the majority 
of SNe IIP, except for rare cases of SN II intermediate between IIP and IIL (linear), of which SN 1970G is an 
example. 

Key words. (Stars:) supernovae: type IIP; Nucleosynthesis, abundances; Photometry, spectra. 



1. Introduction 

Supernovae type II with a plateau in the light curve 
^ (SNe IIP) are thought to be a spectroscopically and 
k> photometrically homogeneous family of core collapse su- 
'vh ■ pernovae HFilippenkolOOTl ). The light curve of SN IIP 
5^ i is characterized by two distinct phases; a plateau with 
a duration of 60 — 100 days and a quasi-exponential 
tail at the later epochs. The plateau corresponds to 
the radiative cooling of the hot opaque envelope orig- 
inally heated by the explosion ( |Grassberg et al. 1971| ). 
Variations of observational properties at the plateau 
phase are interpreted in terms of variation of ejecta 
mass, explosion energy, and radius of the progenitor 
( |Litvinova fc Nadyozhin 1985[ ) . The exponential tail is at- 
tributed to the instant reprocessing of the energy of ra- 
dioactive decay of ^^Co ( [Weaver fc Woosley 1980| ) and the 
luminosity at this stage is directly determined by the mass 
of ejected ^^Ni, a fact successfully exploited for SN 1987A 
dUatchpole et al. T988| ). 

The SN IIP phenomena are believed to originate from 
the explosion of red supergiant stars whose main sequence 
masses lie in the range 10 — 25 Mq; precise boundary 
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are still a subject of controversy, as well as the mech- 
anism(s) of explosion. The study of optical properties 
of SNe IIP provides a means of recovering their major 
characteristics (mass, energy, ^^Ni mass, asymmetry) and 
eventually to impose constraints on the explosion mod- 
els and pre-supernova parameters. The ^^Ni mass is one 
of the crucial parameters since it presumably depends 
on the presupernova structure and the explosion model 
l|Aufderheide et al. 1991|l and can be directly measured. 

The straightforward way to measure ^^Ni mass in 
SNe IIP is based upon the physical argument that the 
early (120-300 d) bolomctric luminosity on the radioac- 
tive tail should be equal to the luminosity of radioac- 
tive decay of ^^Co. This is because at this stage the 
loss of the internal energy to expansion {pdV work) is 
small, while the envelope is nearly opaque to gamma- 
rays. This method was successfully used for SN 1987A 
dCatchpole et al. 1988| [Bouchet fc Danziger 1993| ) with 
estimates 0.08 Mq and 0.07 Mq, accordingly, and an av- 
erage 0.075 Mq. The application of this method to other 
SNe IIP is, however, hampered by the lack of infrared 
observations longward of the I band. Note, the UBVRI 
data account for about half of the luminosity at this age 
according to results for SN 1987A dCatchpole et al. 1988) 
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ISchmidt et al. 1993|l . Therefore, for SNe IIP one uses the 
absolute flux in one band (e.g. V) or several bands in 
the optical, which, on being compared with the absolute 
flux of SN 1987A provides an estimate of ^^Ni mass. This 
method essentially assumes that the absolute flux in one or 
several optical bands is a constant fraction of the bolomet- 
ric flux, i.e. the spectral energy distribution (SED) during 
the radioactive tail epoch is similar for all SNe IIP. Using 
this approach, Phillips et al. (1990) claimed that similar B 
and V absolute magnitudes of SN 1987A, SN 1969L and 
SN 1983K on the radioactive tail implies similarity of their 
ejected ^^Ni masses. A similar conclusion was made by 
Patat et al. (1994) from a study of a large sample of type 
II SNe. Schmidt et al. (1993) shghtly modified this method 
applying a procedure of reconstruction of the light curve 
using a bolometric correction determined from SN 1987A 
data. Using a constant bolometric correction does not nec- 
essarily improve the accuracy of the determination of ^^Ni 
mass compared with the use of absolute magnitudes in 
one or several filters. Recently Hamuy (2003) using this 
method recovered ^^Ni mass for an extensive sample of 
SNe IIP and pointed to an interesting correlation of ^^Ni 
mass with the luminosity (My) on the plateau. 
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Here we revisit the problem of ^^Ni mass in SNe IIP, 
pursuing two major goals. First, we independently deter- 
mine the ^^Ni mass for a different sample using a unified 
approach to distance determination and present a new cor- 
relation between ^^Ni mass and light curve shape. The 
second major goal is to check the possibility of the use 
of Ha luminosity as a tracer of ^^Ni mass. This possibil- 
ity is intriguing bearing in mind that in the case of high 
redshift SN IIP {z > 0.05) we face a problem of contri- 
bution of the host galaxy background, which will inhibit 
the confident estimation of the broad band magnitude at 
the radioactive tail phase. On the other hand, Ha is not 
affected by the stellar background and thus may be used 
as an indicator of the ejected amount of ^^Ni. Nevertheless 
incorrect substraction of any underlying H II region emis- 
sion could hamper accuracy; identification of this narrow 
component in the blend would remove this problem. The 
use of Ha as an indicator of ^^Ni mass is prompted by 
early radioactive models which show an expected propor- 
tionality between Ha luminosity and the amount of ^^Ni 
( |Chugai 1990| ). This has also been noted in connection 
with the comparison of light curve and Ha luminosities 
for SN 1999em and SN 1987A HElmhamdi et al. 2003|l . 

The paper is organized as follows. We first select a 
sample of well observed SNe IIP with late photometry and 
spectra, and adopt appropriate distances and extinction 
(Sect. 2). We then apply the V light curve technique to 
estimate the ^^Ni mass for these supernovae and study 
correlations between the ^^Ni mass and light curve shape 
(Sect. 3). Then we use the late Ha luminosity to derive 
^^Ni mass (Sect. 4). Both photometric and specroscopic 
^^Ni mass are then compared and the implications are 
discussed (Sect. 5). 



Fig. 1. The B — V (upper panel) and V — R (lower panel) 
intrinsic colour evolution of the SNe sample. Both options 
for SN 1997D are shown (sect. 4.2). The used reddening 
is reported in Table 1 (fourth column). 



2. Sample, distance, extinction 

Our basic sample consists of 9 SNe IIP, namely: SN 1987A, 
SN 1969L, SN 1988A, SN 1988H, SN 1990E, SN 1991G, SN 
1992H, SN 1997D and SN 1999em, although throughout 
the paper we will be invoking other objects (SN 1970G, SN 
1995ad, SN 1995V, SN 1995W, SN 1999gi and SN 1999eu). 
The events are selected on the basis of published and un- 
published photometric and spectroscopic observations at 
the plateau and later epochs. The sample parameters are 
presented in Table 1 together with the corresponding ref- 
erences to the literature from which the photometry has 
been taken. In what follows we adopt the standard redden- 
ing laws of CardeUi et al. (1989). For the problem of ^^Ni 
determination the total extinction and distance estimates 
are of principal importance. Galactic extinction is removed 
using the map of galactic dust extinction by Schlegel et 
al (1998). The host galaxy reddening is then estimated 
from the B — V and V — R colour excess compared to the 
intrinsic colour curves of SN 1987A. This approach is jus- 
tified by the fact that at the late photospheric phase SNe 
IIP seem to follow colour evolution similar to SN 1987A 
HSchmidt et al. 1992|l . The use of the B - F colour alone 
to estimate the host galaxy reddening may present some 
problems, as in most cases it provides negative reddening. 
This point was noted also by Hamuy (2003) who used the 
V — I index for an independent estimate of extinction. 
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SN 


P&r6iit 


DistcincG 




A/\/ l&tG slopG 




name 








mag (lOOd)"^ 




1987A 


LMC 


0.05 


0.6 


0.976 


1 2 


1 QQQpm 


NOC 1637 


8.8 


0.31 


0.97 


3 4 5 


1969L 


NOC 1058 


9.057 


0.203 


0.97 


6 7 


1988A 


NGC 4579 


22.95 


0.136 


1.12 


8, 9 


1 yoon 


IN VI V OO 1 O 


9S 4-7 








1990E 


NGC 1035 


16.18 


1.2 


0.86 


11, 14 


1991G 


NGC 4088 


14.11 


0.065 


1.11 


10 


1992H 


NGC 5377 


29.07 


0.054 


0.97 


15 


1997D 


NGC 1536 


16.84 


0.07 


0.87 


12, 13, 16 


1970G 


NGC 5457 


7.2 


0.4 


1.04 


6, 17 


1999gi 


NGC 3184 


10.91 


0.65 




18 


1999eu 


NGC 1097 








19 



REF^ 

1- Arnctt 1996; 2- Hirata et al. 1987; 3- Baron ot al. 2000; 4- Elmhamdi ot al. 2003; 
5- Hamuy ot al. 2001; 6- Kirshner & Kwan 1975; 7- Patat ot al. 1994; 8- Ruiz-lapuonto ot al. 1990; 
9- Turatto ot al. 1993; 10- Blanton ot al. 1995; 11- Bcnotti ot al. 1994; 12- Turatto ot al. 1998; 
13- Bonotti ot al. 2001; 14- Schmidt et al. 1993; 15- Clocchiatti et al. 1996; 16- Zampiori ot al. 2002 
17- Barbon ot al. 1973; 18- Leonard et al. 2002; 19- Pastorello et al. 2003 
* Hn— 70 km s^"^Mpo^^ is adopted. 

The recovered total visual reddening (Ay*) of our 
sample is listed in Table 1. Fig. 1, on the other hand, 
displays the intrinsic B — V and V — R colour evolu- 
tion. Note that for the case of SN 1997D, both colours 
show a large excess by the end of the photospheric phase. 
The small (case b; ITuratto et al. 1998|l and large (case 
a; [Zampieri e al. 2002| ) age scenarios are displayed. This 
high reddening seems unlikely since the presence of in- 
terstellar lines was not reported in the early spectra 
as one may expect in a case with such high reddening 
IjTuratto et a l. 1998). The very red colour at the end of 
the photospheric phase is seen in other objects with a 
very low ejected ^^Ni mass (e.g. SN 1999eu; Pastorello et 
al. 2003, in preparation). An explanation of this peculiar- 
ity may be related to the fast cooling from the plateau 
phase to reach the faint radioactive tail, and/or possibly 
relates to the nature and structure of the progenitor star. 
Interestingly, the late nebular intrinsic B — V colours of 
the complete sample seem to be very similar to that of 
SN 1987A, both in slope and magnitude. The convergence 
to the same fate at later phases is also confirmed for the 
faint object SN 1997D and as well as for the underlumi- 
nous SN 1999eu (Pastorello et al. 2003). Ahhough study- 
ing the spectral energy distribution (SED) of SNe IIP is 
beyond the scope of the present work, the fact of a small 
scatter (~ ±0.3 mag) in the sample around the SN 1987A 
intrinsic {B — V)o colour indicates the similarity in the 
SED for SNe IIP on the radioactive tail. The distance of 
D = 50 kpc is used for SN 1987A, while for the other SNe 
we adopt distances derived from the recession velocity of 
the host galaxy corrected for Local Group infall onto the 
Virgo Cluster as reported in the "LEDA^" extragalactic 

^ available online at ,http: / /leda.univ-lyonl .fr /cgi-bin /single. pl_ accurate value of the distance to SN 1999em, 7.8 Mpc, 
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Fig. 2. The absolute F-magnitude evolution of the SNe 
sample. The used parameters are reported in Table 1. Both 
possibilities for SN 1997D are plotted. The dashed line 
shows the slope for ^^Co decay. 



database. The Hubble constant Ho=70 km s^^Mpc^^ is 
adopted. The computed distances are reported in third 
column of Table 1. We note here that a possibly more 
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Table 2. '^°Ni mass estimates from photometry 



SN 


Mm{V) 


ti 


name 




fdavsl 


1987A 


0.075 


107 (±2) 


1969L 


0.067 (±0.002) 


110 (±4) 


1988A 


0.088 (±0.003) 


134 (±4) 


1988H 


0.033 (±0.004) 




1990E 


0.043 (±0.0024) 




1991G 


0.021 (±0.0032) 


122 (±5) 


1992H 


0.123 (±0.002) 


112 (±5) 


1997D(a) 


0.0065 (±0.0003) 


112 (±5) 


1997D(b) 


0.0036 (±0.0002) 


112 (±5) 


1999em 


0.027 (±0.002) 


112 (±4) 


1970G 


0.051 (±0.003) 


94 (±4) 


1999gi 


0.0246 (±0.0007) 


120 (±3) 


1999eu 


0.0028* 


110 (±4) 


* From PastorcUo ot al. 2003 (In 


preparation) 



comes from stellar studies of the parent galaxy NGC 1637 
HSohn &: Davidge 1998| ). 



3. ^^Ni mass from V light curve 

The luminosity of SN IIP on the radioactive tail is con- 
trolled by the radioactive decay (^^Co— >^^Fe) and if the 
trapping of gamma rays is efRcient, the decline rates of the 
light curves of SNe II should converge to the exponential 
life-time of ^^Co (i.e. 111.26 days or 0.976 mag per 100 
d). The radioactive tails in SN IIP all show similar decay 
rates, especially in the V band (,Patat el al. 1994|l . This 
fact provides a solid basis for using the V light curve for 
the recovery of the ejected ^^Ni mass (Phillips et al. 1990; 
Schmidt et al. 1993). In Fig. 2 we display the absolute V 
light curves of the SNe sample together with that of SN 
1987A. The computed late time decline slopes (150 — 400 
days since explosion) are reported as well in Table 1 (fifth 
column). Their mean value is < 7^ > ~ 0.99 (±0.13), 
consistent with the radioactive decay of ^^Co and conse- 
quent trapping of the gamma-rays. 

The ^^Ni mass is estimated from absolute My mag- 
nitudes between 120 - 400 d using the SN 1987A tail as 
a template. Results produced by the least squares fit are 
reported in Table 2. Both possible cases for SN 1997D are 
reported (more details concerning this event are discussed 
in sect. 4.2). These values show a significant range of the 
ejected ^''Ni masses, from 0.0028 Mq and 0.0036 Mq for 
subluminous SN 1999eu and SN 1997D(case b) respec- 
tively to 0.123 Mq for SN 1992H with the average value 
« 0.05 Mq. Our values differ somewhat from those deter- 
mined recently by Hamuy (2003) for objects in common. 
The differences stem primarily from the different distances 
adopted here. As noted earlier the methodology followed 
by Hamuy is based upon the conversion of V magnitude 
into a bolometric luminosity through an assumed constant 



bolometric correction, while we are adopting SN 1987A lu- 
minosities as a template at late epochs. 

With our values of the ^^Ni masses it is instructive 
to check that the correlation between plateau My magni- 
tude and the amount of ^^Ni found by Hamuy (2003) is 
preserved. We have modified the definition of the plateau 
My magnitude compared to Hamuy (2003) in order to 
avoid any uncertainty in the explosion time. Therefore, 
we use the inflection time during transition from plateau 
to the radioactive tail as a zero point. The inflection time 
ti is defined as the moment when the first derivative at the 
transition phase S = —dMy/dt (we dub it "steepness") is 
maximal. To calculate S we use the following procedure. 
The V band flux in the transition period from plateau and 
radioactive tail is approximated by a sum of plateau and 
radioactive terms: 



F = A 



1 ± (iAo)« 



Scxp(-Vlll.26), 



(1) 



where A, B, to, p and q are parameters derived by the 
minimization technique in the sensitive interval ti ± 50 
days. The behaviour of S and determination of is demon- 
strated for each SN in Fig. O The derived ti are reported 
in Table 2 (third column) with the corresponding errors 
while the computed S values and their errors are presented 
in Fig. 5. The errors in S and ti are estimated according 
to the best fit (shown in Fig. 3) and as well by introducing 
test points for events with poor data especially at the tran- 
sition phase (from plateau to radioactive decay) . With the 
determined inflection time ti we choose the epoch ti — 35 
days as a reference for the My on the plateau. 

Fig. 01 (upper panel) demonstrates the correlation be- 
tween photometric ^^Ni mass and the absolute magnitude 
My{ti — 35) for the sample of SNe IIP with measurable S 
and ti values. This plot shows small scatter around a linear 
trend and thus supports the correlation found by Hamuy 
(2003). In our case, and assuming case "a" (Zampieri et 
al. 2002) for SN 1997D, the linear correlation is described 
by the equation: 



log M(^^Ni) = -0.438My(ti - 35) - 8.46. 



(2) 



Therefore this result suggests a direct correlation between 
V magnitudes (absolute or apparent) of plateau and ra- 
dioactive tail. In fact Fig. 01 (lower panel) demonstrates 
the correlation between My magnitudes at plateau (on 
ti — 35 day) and tail (on ti + 35 day). This is therefore a 
direct form of the correlation shown in the upper panel. 

The sample we explore reveals another interesting cor- 
relation. Fig. |S1 demonstrates that the steepness S anti- 
correlates with ^®Ni mass: the lower is the ^^Ni mass the 
larger is S, i.e. the steeper is the transition from plateau to 
the tail. We note however, that an accurate determination 
of S necessitates a V light curve with a reasonably high 
density of observational points at the end of the plateau 
phase and the beginning of the radioactive tail. Indeed 
SNe 1990E and 1988H are not studied here because of 




Fig. 3. Determination of the steepness and moment of inflection. For each SN of the sample: upper panel displays 
the V light curve (dotted points) together with the best flt (dashed curve) while lower panel shows the steepness S 
(see text). The inflection point ti shown by dashed line corresponds to the maximum of S. SN 1999eu is not shown 
because the data is not yet published. 
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Fig. 4. The correlation between My magnitudes of 
plateau and of radioactive tail. Upper panel shows the 
correlation between My at the moment (ij — 35d) and 
^^Ni mass derived from the tail magnitudes; case "a" of 
SN 1997D is adopted for the fit. Lower panel shows di- 
rectly the correlation of plateau magnitude Mvih — 35d) 
and tail magnitude MviU + 35d). 

the paucity of their light curve data. The best linear fit 
(Fig.[Sl) reads: 

log M(^^Ni) = -6.2295 5 - 0.8147 (3) 

The interpretation of this correlation requires hydrody- 
namical modeling with different amounts of ^®Ni and de- 
grees of mixing. In the case of unmixed ^^Ni one does 
not expect a notable dependence of the steepness on the 
amount of '^^Ni. This is readily seen from the modeling 
SNe IIP hght curves by Eastman et al. (1994). For the 
same amount of ^^Ni the mixing results in a decreased 
steepness because of an increase of radiative diffusion 
IjEastman et al. 1994|) . Unfortunately, in the cited paper 
the authors did not model cases of variable ^^Ni mass with 
similar degrees of mixing, so the question of the physics 
behind the above correlation, remains open. At this stage, 
we may only suggest that somehow the increase of the ^^Ni 
mass in SNe IIP ejecta favours the larger contribution of 
radiative diffusion at the end of the plateau and, therefore, 
a less steep transition from plateau to the radioactive tail. 
It may well be that the increase of the ^^Ni mass is ac- 
companied by the growth of the degree of mixing degree 
which favours a less steep decline as demonstrated in Fig. 
10 and Fig. 12 of Eastman et al. (1994). 



92H 



-1 - ^ 88A 

t*^ 87A 




99eu 



0.1 0.2 0.3 

S (mag/day) 

Fig. 5. The correlation between ^^Ni mass and steepness 
S. Case "a" of SN 1997D is adopted for the linear fit 
(dashed line). 

Interestingly, if the correlation between steepness and 
^^Ni mass is confirmed, this will provide us with an excit- 
ing possibility of probing ^^Ni mass in SNe IIP from the 
light curve shape in a manner independent of the distance 
and extinction. 

4. ^^Ni mass from Ha luminosity 

The Ha luminosity of SNe IIP at late epochs reflects 
the total ionization in the weakly ionized gas caused by 
radioactive decay of ^^Co and, therefore, to a first ap- 
proximation is proportional to the overall deposition rate 
and hence to the ^^Ni ma ss ( |Chugai 1990] IXu et al. 19921 
IKozma Fransson 1992|l unless effects due to higher den- 
sities become important. This suggests that the Ha lumi- 
nosity resulting from gamma-ray trapping may provide a 
quantitative measure of the ^^Ni mass in SNe IIP. The ac- 
curate measurements on SN 1987A provide a control. To 
explore deviations from the simple proportionality we use 
an upgraded radioactive model of the Ha luminosity. 

4.1. Model of Ha luminosity 

Some brief comments about the modified model of Ha lu- 
minosity powered by the radioactive decay are in order. 
The primary purpose of the modification of the previous 
version (Chugai 1990) is to specify better the early nebular 
phase. This was done by the implementation of macro- 
scopic mixing, coUisional and radiative de-excitation ef- 
fects. The supernova envelope is mimicked by two zones. 
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Table 3. Parameters of the Ha evolutionary models 
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Parameter 


Unit 


Standard 


a 


al 


b 


bl 


b2 


c 


d 


e 


f 


fl 


M 




14 


10 


7 


14 


14 


14 


14 


14 


14 


14 


14 


E 




1 


1 


1 


1 


1 


1 


1.5 


1 


1 


1 


1 


Mni 


Mq 


0.075 


0.075 


0.075 


0.0375 


0.00375 


0.15 


0.075 


0.075 


0.075 


0.075 


0.075 


,/mix 




0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.1 


0.4 


0.4 


./"met 




0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.3 


0.6 


Tc 


K 


7000 


7000 


7000 


7000 


7000 


7000 


7000 


5000 


7000 


7000 


7000 



The inner zone with the radius rmix — Wmix^ is a mixing 
core with a total mass fraction /mix, which is divided be- 
tween metals together with He-rich matter (we dub this 
mixture "metals" with the mass fraction of /,n) and H-rich 
matter with the mass fraction 1— /mix- The outer zone con- 
sists of hydrogen-rich matter. All the ^^Ni resides in the 
mixing core and has a generally clumpy distribution. This 
is mimicked by ^®Ni clumps in cocoons of metals with the 
total area of cocoons being 47rr^j;^/s. Here /s is a mixing 
parameter, which is unity if all metals reside in the spher- 
ical layer with the inner spherical cavity occupied by ^®Ni. 
We adopt fg — 2 which corresponds to moderate mixing. 
The escape probability for gamma-rays through the metal 
layer is exp(— T,n), where T,n is the optical depth of the 
metal layer, while the absorption probability in the mix- 
ing zone is described by the expression ti/(1 -f ti), where 
Ti is the optical depth of the inner zone. In the mixing core 
gamma-rays have a finite probability of being absorbed by 
H-rich material, which is proportional to 1 — /mix- The ab- 
sorption probability for the outer hydrogen layer with the 
optical depth T2 is defined as p2 = 1 — exp(— r2). 

The energy deposited in the hydrogen is shared be- 
tween ionization, excitation and thermal energy with 
corresponding branching ratios l|Kozma fc Fransson 19921 
IXu et al. 1992|l . The ionization of hydrogen is calculated 
in the approximation of three levels plus continuum. The 
photoionization from the second level by two-photon and 
continuum radiation was calculated assuming the contin- 
uum luminosity is equal to the deposited luminosity while 
the spectrum is assumed to be dilute black-body with 
the temperature Tc. The electron temperature was set to 
be 5000 K. Effects of radiation transfer in Balmer and 
Paschen continua are treated in the escape probability ap- 
proximation (Chugai 1990). 

The sensitivity of the radioactive model to parameter 
variations (Table 3) is demonstrated in Fig.El Each panel 
shows the effect of a single parameter variation compared 
to a standard model, which is in fact the optimal model 
for the Ha luminosity of SN 1987A. The saturation at 
the early epoch is a new feature compared to the previous 
model and this is related to several factors: gamma-ray 
trapping in the inner zone, collisional de-excitation of the 
third level and Ha absorption in the Paschen continuum, 
especially in the mixing zone. In Fig.|^ we present models 
with lower mass (10 and 7 Mq) compared to the standard 




Days 

Fig. 6. The Ha luminosity in the radioactive model. The 
thick line shows the standard model, while thin line curves 
show models a, b, c, d, e, f. Dotted lines in panels a, b 
and f correspond to models al, bl and fl respectively. 
Additional model is shown in panel b (model b2, dashed 
line). Details of model parameters are presented in Table 
3. 



14 Mq. a lower mass results in a high rate of the luminos- 
ity decay owing to a lower optical depth for gamma-rays. 
For moderate mass variation (factor < 1.4) the Ha lu- 
minosity in the range 200 — 300 d shows less than 10% 
deviation from the standard case. With the rest of param- 
eters fixed the Ha luminosity is proportional to the ^^Ni 
mass (Fig. 6b) at the nebular epoch t > 200 days. This 
property reflects the simple truth that to a good approx- 
imation the total ionization rate of hydrogen in the enve- 
lope is proportional to the total radiactive decay rate, i.e. 
^^Ni mass. The proportionality, however, breaks down for 
t < 200 days because of saturation effects related to Ha 
collisional de-excitation and continuum absorption, which 



8 



Elmhamdi et al.: M{^^m) and Ha in SNe IIP 



become important at early epoch. The effect of factor 1.5 
higher kinetic energy is similar to the effect of factor ^1.2 
lower mass (Fig. IS];), in accordance with the dependence 
of optical depth on both parameters oc M"^ /E. A lower 
continuum temperature (Fig. EJi) results in a lower lu- 
minosity, especially, at the early epoch (t < 200 d) since 
the ionization (and therefore recombination) rate becomes 
lower. A lower mixing core fraction (Fig. results in 
the higher Ha luminosity since the total deposition in the 
outer H-rich matter increases. The higher metal fraction 
in the mixing core (Fig. IHt) , on the contrary, suppresses 
the late time Ha luminosity consistent with the decreased 
deposition into mixed hydrogen material. However, this 
effect is small if the variation of this parameter is within 
a factor 1.5. 

The modeling shows that if ejecta parameters vary less 
than a factor 1.4 compared to those of SN 1987A, the Ha 
luminosity is then proportional to the ^^Ni mass (to within 
10% accuracy) at the late time nebular phase 200 — 400 
days. In this age range using the SN 1987A model as a 
template we may estimate the ^^Ni mass in other SNe IIP 
from the ratio of Ha luminosity to that of SN 1987A. It is 
worth noting here that despite the large possible variation 
of the main sequence mass for SNe IIP, presumably say 
10 — 25 Mq, the range in mass of the ejecta may actually 
be smaller since mass loss increases rapidly with the mass. 

4.2. Results of^^Ni mass determination 

The spectra utilized in this section are based on data in 
the Asiago/ESO SN Catalogue. Additional spectra were 
kindly provided by R. Stathakis. A sample of the spec- 
tra is shown in Fig. 7 together with some line identifica- 
tions. The available spectra of our sample were corrected 
for reddening effects and the recovered integrated Ha line 
fluxes are then translated to luminosities using adopted 
distances. In a straightforward approach we use the Ha 
luminosity of SN 1987A as a template to recover the ^^Ni 
mass from the Ha luminosity of SNe IIP assuming propor- 
tionality between the ^^Ni mass and the Ha luminosity in 
the range between 200 — 400 days. The resulting values of 
^^Ni mass are given in Table 4 in the second column. 

Another approach has been to fit the Ha light curves 
with the model of luminosity evolution. All the other pa- 
rameters, except for ^^Ni mass, are assumed to be the 
same as in the standard model for SN 1987A unless ad- 
ditional information indicates essentially different param- 
eters. The plateau of SN 1970G is very brief [tp « 60 d) 
which indicates a lower ejecta mass than that of a typical 
SN IIP. We adopt 8 Mq as the most appropriate value for 
SN 1970G ejecta taking into account both the light curve 
and the rate of the Ha flux decay. The results are plot- 
ted for 9 SNe in Fig. |S1 Recovered values of ^^Ni mass are 
given in Table 4 (third column) . Note that the straightfor- 
ward approach based upon the use of the Ha light curve 
of SN 1987A (second column of Table 4) and the appli- 
cation of model fitting produce quite similar values with 




wavelength (A) 

Fig. 7. Sample of late time spectra of SNe IIP with the 
corresponding time since explosion. The spectra have been 
corrected for the recession velocities of their host galaxies. 
Some line identifications are also indicated. 



a maximal difference, apart from SN 1970G(SNP/L) and 
SN 1997D(treated bellow), of -13% for SN 1988H. 

This deviation may be adopted as an uncertainty of the 
56]\fj mass determination from the accurate data. Given 
the uncertainty of the observed value of the Ha luminosity 
(at least 10%) we are therefore able to derive a ^^Ni mass 
from Ha with an uncertainty of about 25%. 

The subluminous SN 1997D is a special case. 
Originally, the age of the supernova at discovery has been 
estimated to be — 50 days (Turatto et al. 1998), while 
Zampieri et al. (2002) argue for an age almost twice as 
large ~ 95 days. We consider both large and small 
age options. In the large age case we adopt M ~ 18 Mq, 
E = 9 X 10^° erg (Zampieri et al. 2002), while in the 
small age option we adopt M = 6 M0 and E = 1 x 10^" 
erg ( |Chugai &: Utrobin 2000| ). A reasonable fit in the first 
case (Fig. Et) is found for the ^^Ni mass of 0.011 Mq 
with a mixed core fraction /mix = 0.05, which corre- 
sponds to the mixing core having 0.9 Mq and a core 
metal fraction /m = 0.99, i.e. with the metal core practi- 
cally devoid of hydrogen. The model also requires a low 
continuum temperature in the ultraviolet Tc = 4500 K. 
For the small age choice the fit (Fig. Eb) corresponds 
to M(56Ni) = 0.0024 Mq. In this case /^ix = 0.05, 
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Fig. 8. The Ha luminosity evolution in different SN IIP. 
Models {solid lines) are overplotted on the observational 
data dots (sect. 4.2). 

/m = 0.3 (a significant amount of mixed hydrogen) and 
Tc = 5500 K. 

The ^^Ni mass derived from the Ha luminosity adopt- 
ing SN 1987A as template, MNi(Ha), is plotted in Fig.lTIil 
versus photometric '^^Ni mass, M-^i{V), derived from the 
tail Mv magnitude (Table 2). Both sets of values agree 
within 20% . This consistency supports the proposition 
that the Ha luminosity may be a good indicator of ^^Ni 
mass in SNe IIP unless we are dealing with extreme cases 
such as SN 1970G (SNP/L) or underluminous cases such 
as SN 1997D. 

4.3. Application of Ha to ^^Ni diagnostics 

In the absence of the late time photometry but with avail- 
able spectra of SNe IIP at the nebular epoch we may use 
the Ha luminosity to estimate ^^Ni mass ejected by super- 
nova provided the explosion time is known. We demon- 
strate this approach to SN 1995ad, SN 1995V and SN 
1995W (all type IIP) for which late photometric data are 
not available. Distances are computed using the corrected 
recession velocities (reported in "LED A" data base assum- 
ing Ho =70 km s-^Mpc"^). 

SN 1995ad: This SN was discovered on Sep. 28.8 UT 
in NGC 2139 by R. Evans IjlAUC 623911 . Its Sep. 29.3 UT 
spectrum, obtained by S. Benetti, displays a hot contin- 
uum with temperature Tbf, ^ 13000 K and a broad Ha 
emission {FWHM =11000 km s^^). These characteris- 
tics combined with the fact that nothing was seen at the 



Table 4. Recovered mass of ^^Ni using SN 1987A as 
template (second column), and from late time L(Ha) 
modeling "M^(iJa)". 



blN 




71 j-m /TT _ \ 


name 






1969L 


0.068(±0.004) 


0.065 


1970G 


0.029(±0.011) 


0.035 


1988A 


0.085(±0.004) 


0.083 


1988H 


0.039(±0.005) 


0.045 


1990E 


0.04(±0.003) 


0.038 


1991G 


0.022(±0.0005) 


0.025 


1992H 


0.106(±0.008) 


0.105 


1997D(a) 


0.0028(±0.0006) 


0.011 


1997D(b) 


0.0018(±0.0004) 


0.0024 


1999em 


0.022(±0.003) 


0.024 



position of the SN in the Aug. 25 image (lAUC 6239) 
provide constraints on the explosion time. The assumed 
distance and reddening are respectively D — 23.52 Mpc 
and Av = 0.112 mag (NED; Schlegel et al. 1998) while 
the galaxy inclination is 40.9°. The spectrum we use here 
is taken 1996 March. 24, about 180 days since explosion. 
The derived Ha luminosity, ^ 1.06 x 10^° erg s~^, is then 
compared to that of SN 1987A at a similar epoch to re- 
cover an ejected mass « 0.056 AIq. 

SN 1995V: discovered in NGC 1087 by R. Evans 
flAU C; 619711 . The galaxy inclination is 33.2°. The 
adopted distance and reddening are, accordingly, D — 
20.61 Mpc and Ay = 1.37 mag, while the explosion 
time is assumed to be 1995 July. 25 IjFassia et al. 1998jl . 
The available late spectrum was taken around day 409 
since explosion. Comparing the recovered Ha luminosity, 
- 1.33 X 10^^ erg s-\ with that of SN 1987A leads to the 
estimate of Mi^^Ni) « 0.046 Mq. 

SN 1995W: discovered on Aug. 5.65 UT in NGC 
7650 by A. WiUiams and R. Martin ( HAUC 6206|l . The 
spectrum taken on Aug 17.29 UT displayed features of 
SN IIP around one month (lAUC 6206). The distance is 
D — 44.81 Mpc, and the reddening is Ay — 0.044 mag 
(NED; Schlegel et al. 1998) and the host galaxy inclina- 
tion is 47.4°. The late spectrum was taken around day 
300 (1996 May. 12). The Ha luminosity, - 4.79 x lO^^ erg 
s~^, when compared to that of SN 1987A at similar epoch 
yields an amount of ^^Ni mass of M{'^^lS!i) « 0.048 M©. 
Note that the adopted reddening values for SN 1995ad 
and SN 1995W do not include the host galaxy reddening. 
For these events the inclination of their host galaxies may 
provide further indication of uncertainty in the extinction 
since large extinctions might be expected for highly in- 
clined galaxies if the SN occurred towards the far side. 
Although no sign of high extinction from interstellar lines 
was reported in early spectra of these two objects, the de- 
rived amounts of Af (^^Ni) might be taken as lower limits 
bearing in mind the significant inclination of their host 
galaxies. 
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Fig. 9. Modeling Ha luminosity for SN 1997D. The up- 
per panel shows the large age option, while lower panel 
shows the small age option. Model results (solid lines) are 
overplotted on observational dots. 



5. Discussion and conclusion 

The primary goal of this paper was to demonstrate that 
the Ha luminosity at the nebular epoch may be an indi- 
cator of the ^^Ni mass in SNe HP ejecta. To explore this 
idea we selected a sample of well observed SNe HP, and 
using V magnitudes on the radioactive tail we derived 
in a standard way the photometric mass of ^^Ni mass. 
With our sample we confirmed using a slightly different 
approach the correlation between ^^Ni mass and absolute 
magnitude My of the plateau reported recently by Hamuy 
(2003). We then applied a two-zone model of the Ha lu- 
minosity in SN HP to explore the sensitivity of the Ha 
behaviour to variation of model parameters. We found 
that if mass, energy and mixing conditions do not vary 
strongly among SNe HP (less than factor 1.4) then with an 
accuracy better than 10% Ha luminosity is proportional 
to ^®Ni mass during the 200 — 400 days after explosion. 
Ha luminosities were then used to derive ^^Ni masses. 
This was done employing two approaches: first, using the 
Ha light curve in SN 1987A as template and, second, ap- 
plying the model computations. Both approaches agree 
within 15% unless we are dealing with extreme cases such 
as SN 1970G (type IIP/L) and underluminous SN 1997D. 
In both these cases we should possess additional informa- 
tion about ejecta mass and energy to derive the '"'^Ni mass 
from Ha. The ^^Ni mass values derived from the photome- 
try and Ha luminosity agree within 20%, which thus gives 
us confidence that Ha is a good indicator of the amount of 
^^Ni in SNe HP. Simultaneously, this consistency suggests 
that parameters of SNe HP (mass, energy and mixing) 
are not very different. In fact this is consistent with the 




MNi(V) (MJ 



Fig. 10. The correlation between the ^^Ni mass de- 
rived from the Ha luminosity, MNi(Ha), and that de- 
rived from the tail My magnitude, M^i{V). The dashed 
line has a slope of unity. Clear deviation is seen for SN 
1970G(SNP/L) and for both scenarios of the faint event 
SN 1997D. 



uniformity of plateau luminosities and plateau lengths of 
SNe IIP. 

We applied the method of ^^Ni mass estimation 
from Ha for three SNe IIP (SN 1995ad, SN 1995V and 
SN 1995W) without photometry at the nebular phase and 
derived ^^Ni masses of « 0.056 M©, 0.046 Mq, 0.048 Mq, 
quite reasonable values, although they could be underesti- 
mates, since the host galaxy extinction (for SN 1995ad & 
1995W) has not been taken into account. Nevertheless, 
this is a good demonstration of the possibility of the 
method. Generally, the approach based upon Ha may be 
indispensable in cases, when the photometry at the neb- 
ular epoch is absent, or when there is a problem with 
subtraction of stellar background (SN IIP in the bulge, or 
in high redshift galaxies). 

An interesting by-product of the analysis of the 
SNe IIP sample is the demonstrated correlation between 
^^Ni mass and the steepness parameter (S) introduced to 
measure the light curve decay rate at the inflection point. 
The correlation is such that the steeper the decline at the 
inflection point the lower is the ^^Ni mass. Thus radia- 
tive diffusion times and ^^Ni masses are linked. How an 
increased radioactive energy input leads to a higher effec- 
tive opacity will require elaboration by hydrodynamical 
modeling. This correlation, if confirmed, will provide us 
with a distance and extinction independent way to deter- 
mine the amount of ejected ^^Ni. 
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Finally we note that in Fig. 10 the clustering of the 
points around two values of ^^Ni mass viz. 0.005 and 0.05 
may result from poor statistical sampling. On the 
other hand it may be a hint that a mechanism such as 
fall-back is an important one in the evolution of the low- 
mass group. 
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